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We present a study of resonant vibrational coupling between adsorbates and an elastic substrate at low macroscopic 
coverages. In the first part of the paper we consider the situation when adsorbates form aggregates with high local 
coverage. Based upon our previously published theory, we derive formulas describing the damping rate of adsorbate 
vibrations for two cases of such aggregation: (i) adsorbates attached to step edges and (ii) adsorbates forming two- 
dimensional islands. We have shown that damping is governed by local coverage. Particularly, for a wide range of 
resonant frequencies, the damping rate of adsorbates forming well separated islands is described by the damping rate 
formula for a periodic overlayer with the coverage equal to the local coverage in the island. The second part of the 
paper is devoted to facilitating the evaluation of damping rates for a disordered overlayer. The formula describing the 
damping rate involves the parameter (3 which is related to the local density of phonon states at the substrate surface 
and does not allow a closed-form representation. For substrates of isotropic and cubic symmetries, we have developed 
a good analytical approximation to this parameter. For a vast majority of cubic substrates the difference between the 
analytical approximation and numerical calculation does not exceed 4%. 
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I. INTRODUCTION 

Molecules or atoms adsorbed on surfaces of materials can oscillate about the adsorption bond. If this motion 
corresponds, in the absence of the adsorption bond, to free translation parallel to the surface, then this mode is called 
a frustrated translation (FT). It plays an important role in such dynamical processes as reactivity, surface diffusion, 
and desorption. For many systems, especially with a metallic substrate, the frequency of the FT mode is so low that 
it readily couples to low-energy excitations in the substrate, which in the case of metals include both phonons and 
electron-hole pair excitations. These bulk excitations provide decay channels for FT modes, which thus acquire a finite 
lifetime t and become resonances. Associated with this lifetime is the resonance width, or damping rate, 7 = IjIt. In 
many cases the electronic contribution to the damping rate is small compared with the phononic contribution.^ 

We recently developed a general theory of phonon-mediated FT relaxation that can be applied to any adsorbate 
overlayer pattern. Using this theory we derived simple FT damping rate laws for three important classes of adsorbate 
patterns: an isolated adsorbate, an ordered overlayer, and a randomly disordered overlayer. Although there are no 
experimental results to. compare to for isolated adsorbates, the agreement with experiment for ordered and disordered 
over layers is excellent B 

However, these three structural models do not cover all possible experimental conditions. Depending on the inter- 
actions of the adsorbates with the substrate and with each other, their distribution on the surface can differ strongly 
from the periodic or random uncorrelated cases. Under certain conditions the adsorbates may gather into aggregates 
even at very low coverages. In this paper we consider two possible examples of such adsorbate aggregation: (i) adsor- 
bates accumulating along the edges of surface steps, and (ii) adsorbates forming two-dimensional (2D) islands. We 
assume that the average coverage is low enough that these aggregates are well separated from each other and can be 
treated as isolated. For each of these cases we derive a corresponding damping rate formula from our general result. 

Another goal of this paper is to facilitate the evaluation of damping rates for a disordered overlayer. The main 
difficulty here is that one needs to integrate the Green's function for a semi-infinite clastic medium over the in-plane 
wavevector at a dense array of frequencies. Regrettably, no closed- form expression for this integral exists, and it 
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must be computed numerically. In this paper we present a good analytical approximation of this integral for elastic 
media of isotropic and cubic symmetries. We have tested our approximation formula for 88 different materials with 
cubic symmetry, comparing it against well-converged numerical calculations. The difference between the analytical 
approximation and the reference numerical values never exceeds 6%, and for a vast majority of the materials is under 
4%. 

The remainder of this paper is organized as follows. In Section || we briefly review our theory of phonon- mediated 
FT damping. For details of the derivation, the reader is referred to Ref. g. This theory is then applied in Sections III 
and IV to the cases of adsorbates aggregated along step edges and in islands, respectively. In Section |v| we present 
our approximation formula for the /c-space integral of the substrate Green's function, thus facilitating the evaluation 
of the damping rate of a disordered overlayer. Finally, in Section [v| we summarize and draw conclusions. 



II. REVIEW OF THEORY 



The phonon contribution to the damping of law-frequency adsorbate vibrations was considered earlier by Persson 
and Rybergcl (PR) and by Hall, Mills, and Blacka (HMB). Persson and Ryberg considered a single, isolated adsorbate 
coupled to a semi-infinite isotropic elastic substrate. Since direct adsorbate-adsorbate interactions are weak for all 
but the densest coverages, this case was presumed applicable for most experimental conditions. Their model gives a 
damping rate law that varies with the vibrational frequency as 0Jq. However when applied to various experimental 
systems, the PR model predicts damping rates that are much smaller than the measured values. Furthermore, the PR 
model has no explicit coverage dependence; coverage enters the damping-rate law only implicitly through its effect 
on the resonance frequency. The applicability of the model is limited by extremely low coverage (< 1%). Hall, Mills, 
and Black considered a periodic overlayer of rare-gas molecules physisorbed at an isotropic substrate. The damping 
rate implied by their result is proportional to the coverage and to the square of the resonant frequency. The periodic 
distribution of the adsorbates over the surface implied by the HMB model generally occurs at high coverage (> 25%). 
In general, most adsorbate systems have structures that differ from both of these special cases. 

In our recent workEI we developed a general theory for a completely arbitrary overlayer structure. We considered an 
array of N non-interacting adsorbates of mass m located at arbitrary positions {R a } at the surface of an anisotropic, 
semi- infinite elastic substrate of mass density p and elastic modulus tensor Cijki, occupying the half-space z < 0. The 
adsorbates are coupled to the substrate by harmonic springs of frequency u> Q and are allowed to oscillate along one of 
the Cartesian axes a. The adsorbates do not interact with each other directly. We showed that the damping rate for 
this system is described by 
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where Dij(k,uj) is the in-plane Fourier transform of the substrate elastic Green's function evaluated at the surface 
(z = z' = 0), and S(k) is the structure factor defined as 



s & = it ( E ex ?H ■ (R a - Rb)} ) , (2) 




where the angle brackets denote the appropriate ensemble average. The integration in Eq. (|l|) is taken over the bare- 
surface Brillouin zone. The integral I{u>) in Eq. ([!]) is a convolution of two quantities, one of which (inside the square 
brackets) describes all the substrate properties, while S{k) provides all the information on how the adsorbates are 
distributed over the surface. Except for a unitless factor of order unity, the quantity in the square brackets is equal 
to the k- resolved density of elastic modes at the substrate surface for frequency ujq. Any possible explicit dependence 
of the damping rate upon the coverage or the structural order comes from the structure factor S(k). 

We previously applied this general result to three different structure models: an isolated adsorbate, a periodic array 
of adsorbates (appropriate at high coverages), and a random array of adsorbates (appropriate at low coverages). For 
each of these models, the structure factor has a simple algebraic form. Their substitution into Eq. ([j]) leads to 
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for an isolated adsorbate, 
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for a periodic overlayer, Orel and 



7 =(l-0)^ Wo 4 /?+-u; o 2 — ' (5) 

for a random overlayer,i where is the coverage, Aq is the area per surface site, is the speed of an a-polarized 
sound wave propagating perpendicular to the surface, and j3 is defined according to 

1 r d 2 k 
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P=- I 7^2 ^D aa (k,u) (6) 



The integral in Eq. (||) is proportional to the density of a-polarized states divided by frequency. For a 3D elastic 
medium the density of a-polarized states is proportional to w 2 , and thus /3 is independent of frequency to. The 
evaluation of this parameter is discussedJn Section |v|. The damping rates given by Eqs. (||) and (|) are just the 
anisotropic versions of the PRE and HMBq results, respectively No experimental results are available for an isolated 
adsorbate. However, we compared the other two models with experiments on CO on Cu(OOl) by applying Eq. (Q) ta 
the case of an ordered overlayer at 50% coverage,!] and Eq. (|5|) to the case of a disordered overlayer at 3% coverage.Q 
Both models predict the lifetimes in excellent agreement with the measured values. Equation rtq) was independently 
derived later by Persson et a/£3 for the case of an isotropic substrate. They also considered the case of an isolated 
adsorbate of finite size. 



III. ADSORBATES ALIGNED ALONG STEP EDGES 



In this Section we study the case of adsorbates located at step edges. This geometry is important since steps are 
always present at surfaces and they generally attract adsorbates to their edges.Eil While the macroscopic coverage is 
low for this case, the local coverage is high. 

In our analysis we ignore the influence of the steps on the substrate Green's function. The only effect of steps that 
we consider is to set the distribution of adsorbates. We assume that the adsorbates form one-dimensional (ID), locally 
periodic arrays. Starting from Eq. ([l]) we derive the damping rate for an infinite, ID periodic chain of adsorbates, 
and then we extend the derivation for a finite, locally periodic chain. 



A. One-dimensional infinite chain 



We assume that the adsorbates are arranged in an infinite, ID, periodic array along the x axis. The adsorbate 
positions are given by R a = (X a , 0), and the nearest-neighbor distance is d. The ID coverage of the step edge can be 
defined as 6\d = cio/d, where ag is the nearest- neighbor distance between the substrate atoms along the step edge. 
Then the structure factor S(k) takes the form 

n 

where {G n } is the set of ID reciprocal vectors defined as G n = 2itn/d. Substitution of Eq. (^) into Eq. (|l|) leads to 

7 = I If £[ Im A* Q (G„,A^o)] • (8) 

P n 

Since there are no resonant phonon states available in the substrate for in-plane wavevectors \G n \ > uq/cr, where cr 
is the speed of the Rayleigh wave, the imaginary part of the Green's function in Eq. (||) vanishes for these values of G n , 
and thus these terms do not contribute io the sum. Similar to the case of a two-dimensional (2D) periodic overlayer 
discussed in detail in our earlier work,Q it turns out that for all ID coverages presently accessible to experiment 

{Qid > §ad ~ 10%), only the term with G n — contributes to the damping rate. Moreover, the integral of this term 
in Eq. (||) does not depend on frequency ujq, and we can define a constant £ Q according to 
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It turns out that Cx (vibration parallel to the step) differs significantly from £ y (vibration perpendicular to the step). 
However, for substrates with isotropic or cubic symmetry, these values, while different from each other, are independent 
of the orientation of the step. It is, therefore, appropriate to use the notation £|| for oscillations along the step and Cj_ 
for oscillations perpendicular to the step but still in the plane. We can now write the damping rate formula compactly 
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Thus, the oscillations along the step are damped at a different rate than the oscillations perpendicular to the step. 
Typically, the damping rate is higher for vibrations along the step. 

Experimentally measured damping rates can be related to Eq. ( flO| ) in two possible ways. If the frequencies for both 
vibrational modes are close, then the measured damping rate will correspond to an average of 7|| and 7j_. Otherwise, 
both modes can be resolved experimentally, each exhibiting its own damping rate. 



B. One-dimensional finite chain 

Here we consider a finite, locally periodic chain of N adsorbates along one of the step edges. As before, the x-axis 
is along the step, and the nearest-neighbor distance in the chain is d. The total length of the chain is L = N d. The 
structure factor corresponding to this configuration is given by 
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Most of the weight of this function is in the peaks centered at the reciprocal vectors G n — 2nn/d. Even at relatively 
small values of N the cross terms in Eq. ( |j"l| ) may be neglected, and the square of the sum can be approximated by 
a sum of squares 



sin [(k x - G n ) L/2] 
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Equations (|llj) and (|12j) are compared for a chain of N = 6 adsorbates separated by d = 5.54 A(twice the nearest- 
neighbor distance for platinum) in Fig. |l|. The solid line corresponds to the structure factor computed according to 
Eq. (|j"l|), while the dashed line represents S(k) given by Eq. (12). As N increases, the main peaks become narrower 
(the width of the peaks is approximately 2w/L), and the difference between the two formulas decreases. In the limit 
N — > co, both expressions take the form of Eq. (^). 

Finally, we show that the damping rate for a finite chain of the adsorbates can be computed approximately from 
Eq. (|l0|) for large enough N (where "large enough" may be as small as 6). First, we assume that all the weight in the 
structure factor comes from the main peaks. Then, repeating the argument we gave for the infinite chain we conclude 
that if Wo/cr < 2tt (1 — l/N) / d, only the peak centered at G = contributes to the integral in Eq. ([!]), and the 
damping rate becomes 
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If the peak is narrow enough, the imaginary part of the Green's function does not change significantly over its width 
(the region k x < 2it/L). Thus, for each value of k y the Green's function can be replaced by its value at k x — 0. The 
integral over k x can now be carried out explicitly, and Eq. (|l^) follows. 



IV. ADSORBATES FORMING ISLANDS 



In this Section we analyze the damping of the FT mode when adsorbates form 2D islands at the substrate surface. 
As with ID step-edge aggregation, this case corresponds to low macroscopic coverage, but high local coverage. We 
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FIG. 1. The structure factor S(k x ) for a chain of 6 atoms separated by d = 5.54 A. The solid (dashed) curve was computed 
according to Eq. (|ll|) (Eq. @). 



assume that each island is a locally periodic array of N adsorbates with local coverage 9\ oc and area A. We also 
assume that the distance between islands is large enough that the islands can be treated as isolated. 

Unlike in the previous Section, there is no need here for considering the case of an infinite island, because an infinite 
island is the same as a periodic overlayer with 9\ oc — 9, and its FT damping rate is given by Eq. (0). The problem 
of a finite island appears to be more complicated than that of a finite ID chain at a step edge. The reason is that, 
apart from its local coverage and size, the island is also characterized by its shape. It is convenient to represent the 
shape of the island by a function f(r) defined as 

_ J 1 if r is inside island , 
J ^ ' \ otherwise ^ ' 

It is straightforward to show that the structure factor can be expressed in terms of the Fourier transform f(k) of this 
function as 

S{k) = -pY,f{k-G)f{k-G'). i (15) 

G,G' 

where summation is over reciprocal vectors defined by the local periodicity of the adsorbate distribution inside the 
island. 

We further assume that the island is a convex set with similar dimensions in all directions (i.e. its width is of the 
same order of magnitude as its length) . If the number of adsorbates in the island is not too small (the meaning of 
which is discussed below), then the function f(k) will have a strong peak at k = 0, with most of the function's weight 
within the peak. Analytical examples of this function include 

t _ sin [k x a/2] sin[fc y 6/2] 
Mk) - A [k x a/2] [k y b/2} ■ (16) 

for a rectangular island with dimensions a along the x-axis and b along the y-axis, and 
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for a circular island of radius R, where J\ (x) is the first-order Bessel function. 

We now follow the arguments of the ID case. First, we neglect the cross-terms in Eq. (|l^) and write it as 



s® = ^E 3 ■ ( 18 ) 
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Then, we assume that all the weight in the structure factor comes from the main peaks. If the resonant frequency of 
the FT mode u>q is small enough that luq/cr < Gi — A , where G\ is the magnitude of the smallest reciprocal vector, 
and A is the maximal width of the main peak in S(k), then only the G = term will contribute to the integral in 
Eq. (|l|). Thus this equation now takes the form 
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If the the peak of [f(k)] 2 is narrow enough, the imaginary part of the Green's function does not change significantly over 
its width (i.e., the region |fc| < A), and it can be replaced by its value at k — given by ImD act (0, u)q) — l/(c " u>o) , 
and taken out of the integral. The remaining integral can be evaluated with the help of the Parseval's theorem as 



(2tt) 2 A 2 

where, as before, Aq is the area of the surface unit cell. Substituting this result into Eq. (|l9|) we obtain 



m 2 <71oc , on x 

7 ^7^' (20) 

which is exactly the damping rate of a periodic overlayer given by Eq. (^) with the total coverage 6 being replaced 
by the local coverage 9\ oc . 

Finally, we need to determine for which values of N the approximation given by Eq. (^o|) really works. We consider, 
as an example, an isolated island of hexagonal shape on the Pt(lll) surface with local coverage (^oc = 1/4. This is 
about the lowest local coverage observed for islands of small molecular adsorbates. If the approximation given by 
Eq. (^Oj) holds for this coverage, it will also hold for higher coverages. We have computed the integral I(u>) in Eq. (|l|) 
for hexagonal islands of different sizes, and, in Fig. ||, plot I{ui) vs. u> for each island size. The value of the integral 
for the two limiting cases of an isolated adsorbate and an infinite periodic overlayer are shown by the dashed and 
solid lines, respectively. The dotted, dash-dotted, and dash-double-dotted lines correspond to the cases of one, two, 
and three hexagonal shells of adsorbates around the center adsorbate, respectively. The inset to Fig. || characterizes 
the islands by their circumscribing radius R, and the number of adsorbates, N, within the island. 

The behavior of I(u>) shown in Fig. ^|can readily be explained in terms of resonant phonon wavelengths at frequency 
lo. For small frequencies the phonon wavelengths are much greater than the island size, and the adsorbates in the 
island behave like a single isolated adsorbate of mass Nm. Thus, for small frequencies we have I(u>) = Nfiux 2 in accord 
with Eq. (§). Mathematically, this result can be obtained from Eq. ( [l9| ) by setting f(k) — /(0) = A. Frequencies 
for which the resonant phonon wavelengths become comparable with the island size correspond to the cross-over 
region where this quadratic behavior changes. When the phonon wavelengths are smaller than the island size, but 
still significantly greater than the nearest-neighbor distance between the adsorbates, only phonons with in-plane 
wavevector component, k, near zero couple strongly to the FT mode due to the G = peak in the structure factor 
given by Eq. (|l8|). Terms with G ^ have no corresponding resonant phonons at these frequencies. In this frequency 
region we expect I(uj) to follow closely its limiting value of a full periodic overlayer. Finally, when to is large enough 
that the phonon wavelengths become comparable to the inter-adsorbate nearest-neighbor distance, non-zero G-vectors 
in the structure factor become available channels for resonant phonon coupling. However, for realistic local coverages 
($ioc > 1/4) these wavevectors lie at the boundary of the bare-surface Brillouin zone, where continuum models are no 
longer applicable. For the case 8\ oc = 1, this second cross-over region does not exist, because all non-zero G-vectors 
lie outside the surface Brillouin zone. 

From Fig. ^| one can determine the range of resonant FT frequencies for which the approximation given by Eq. ( pp| ) 
holds. Even for the island with only one hexagonal shell (N = 7), there is such a frequency range: I(ui) stays close 
to its limiting value of a periodic overlayer for a range of frequencies between about 3 and 7 meV. For larger islands 
this range expands to frequencies between 2 and 11 meV. 

We can analytically estimate this frequency range for an island of size D. The low-frequency threshold corresponds 
to a transverse wave with wavevector magnitude equal to the width, A, of a peak in the structure factor, which can 



G 



0.025 



0.020 



a> 0.015 



0.010 



0.005 



0.000 



R 



R 



= 1/4 




0a (N = 


1) 


2a (N = 


V) 


4a (JV = 


19) 


6a (N = 


37) 


oo (N-> 


») 



/ 

/'/ 

/ ' //I 

i !, 

/ / 

< // 
_ - ' ./ / 



/ 



5 10 

Frequency [meV] 



15 



FIG. 2. Function /(w) for an isolated adsorbate (dashed line), for an infinite periodic overlayer (solid line) at coverage 
6 = 1/4, and for one-shell (dotted line), two-shell (dash-dotted line), and three-shell (dash-double-dotted line) hexagonal 
islands at local coverage 9\ oc — 1/4. The substrate is Pt(lll). In the inset, R is the circumscribing radius of an island, a is the 
nearest-neighbor distance between the substrate atoms, and N is the number of adsorbates inside an island. The features at 
13.3 meV are caused by coupling to Rayleigh waves at Gi . 



be estimated as 2ir/D. The frequency of such a wave is ui\ ow — 2itct/D. The high frequency threshold corresponds 
to a transverse wave with wavevector magnitude equal to G\ — A, which is approximately 2 {D/d — 1) ttc^/D, where 
d is the inter-adsorbate nearest-neighbor distance, as before. The frequency of this wave is w^igh = {D/d — l)o;i ow . 
For a three-shell hexagonal island (N = 37) with 9\ oc = 1/4 on the Pt(lll) surface, we have d = 5.54 A, D/d = 6, 
and ct = 11 A meV, which leads to the estimates lu\ ow = 2.1 meV and Whigh = 10.5 meV. This frequency range is in 
excellent agreement with the corresponding plot in Fig. ^. Thus, the example of a hexagonal island on the Pt(lll) 
surface demonstrates that, when most of the adsorbates form widely separated islands, the damping rate can be 
computed from Eq. (20) for FT resonant frequencies within a broad range. 



V. APPROXIMATION FORMULA FOR /3 PARAMETER 

The main difficulty researchers may encounter in using Eq. ^ to determine the damping rate for low-coverage 
systems is the necessity to compute the substrate Green's function and the integral in Eq. (|fj) numerically. In this 
Section we present an analytical approximation to the parameter (3 defined by Eq. (|^) . The approximation is developed 
for the cases of systems with isotropic or cubic symmetry. 

We start from the observation that [3 must be a function of the substrate parameters of the model, namely, the 
mass density p and the elastic constants Cy. Generally, i,j — 1,...,6, but we include only independent constants 
characteristic of a given symmetry. This function must obey certain global constraints, however. It follows from the 
equation of motion and the boundary conditions that any quantity depending on the model inputs may depend only 
on the ratios Cij/p, which we denote as fey. Because (3 does not depend on the frequency u>, one can show from the 
equations of motion that (3, as a function of density-normalized elastic constants fcy , must satisfy a scaling property 
given by 

(3(vhi) = v-^ 2 f3(k l3 ) , (21) 
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FIG. 3. Fitting function g(x) defined in Eq. Q23) , computed numerically (squares) and using the analytical approximation 
given in Eq. ( p^ ) (solid line). 

where v is an arbitrary number. This scaling property reduces the number of independent variables in the function 
(3(kij) by one. We have found an excellent analytical approximation to this function for the cases of isotropic and 
cubic substrates. 

A. Isotropic substrate 

An isotropic medium is characterized by only two density-normalized elastic constants: k\\ and k^. The transverse 
and the longitudinal speeds of sound are related to_ these constants by ct = \Jk^ and cl = y/ku, respectively. Taking 



(22) 



into account the scaling property given by Eq. (|2lj), we can express /3 as 

r 2 = kl 4 g(^ 
\ fc 44 



where g(x) is an unknown dimensionless function of one variable. 

To elucidate the behavior of g(x), we have computed /3 numerically for a wide range of kn/k^. The result of these 
calculations is shown by squares in Fig. |[ From Fig. || we obtain that (i) g(x) goes to zero at x = 1 as (x — 1)", 
where 1 < a < 2, and (ii) g(x) asymptotically approaches a constant as x — * oo. We find that the function 



g(x) = P 



(x - If' 2 
3/2 + {x- l) 3 / 2 



(23) 



almost perfectly approximates the numerical data (solid line in Fig. |^). The parameter P is dimensionless and has an 
optimal value of P = 9.14 x 10~ u . 



B. Cubic substrate 



The isotropic approximation is not correct for most single-crystal materials. We therefore extend our approximation 
formula for (3 to describe substrate materials with cubic symmetry, since these are the most commonly used in 
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TABLE I. Density-normalized elastic moduli (10 10 cm 2 /s 2 ), a f3 n (computed numerically), /3 a (computed using the analytical 
approximation), and (/3 a — /3 n )//3n for several cubic materials. All values of (3 are given in 1CP 5 meV _3 A -3 . 



Crystal 


fcn 


k\2 


&44 




/3a 


(A - AO /A 


Al 


39.58 


22.40 


10.51 


1.3635 


1.3657 


0.16% 


Cs 


1.25 


1.04 


0.75 


242.75 


231.03 


-4.83% 


Cu 


18.84 


13.67 


8.48 


3.6321 


3.5586 


-2.02% 


GaAs 


22.34 


10.10 


11.18 


1.7750 


1.7835 


0.48% 


Ge 


24.16 


9.08 


12.55 


1.4659 


1.4608 


-0.35% 


Au 


9.98 


8.45 


2.18 


23.549 


22.361 


-5.05% 


Ir 


25.75 


10.75 


11.37 


1.5337 


1.5377 


0.26% 


Fe 


28.73 


17.80 


14.74 


1.4481 


1.4374 


-0.74% 


MgO 


83.01 


26.64 


43.62 


0.21876 


0.22036 


0.73% 


Ni 


27.85 


17.38 


13.94 


1.5574 


1.5439 


-0.87% 


Pd 


18.87 


14.62 


5.96 


5.3473 


5.1660 


-3.39% 


Pt 


16.13 


11.66 


3.56 


7.9305 


7.7329 


-2.49% 


Rb 


1.87 


1.58 


1.08 


142.05 


133.94 


-5.71% 


Si 


71.12 


27.43 


34.16 


0.29820 


0.30818 


3.35% 


Ag 


11.81 


8.92 


4.39 


9.0756 


8.7958 


-3.08% 



a Mass density and elastic constants were taken from Ref. 



experiments. Crystals with cubic symmetry have three independent density-normalized elastic constants: fcn, fci2, 
and &44. The isotropic limit can be obtained from cubic symmetry by setting k\2 = fcn — 2fc44. Taking advantage of 
the scaling property given by Eq. plf), we can write 



fcl2 fcl2 



(24) 



where h(x, y) is an unknown function of two variables. Generally, the form of the function h(x, y) may depend on the 
particular choice of the surface indices. However, numerical calculations of (3 for different low-index surfaces of the 
same material show that f3 is essentially independent of the surface orientation with respect to the crystallographic 
axes. Thus, the same function h{x 1 y) works well for any low-index surface. 

First, we study how h(x 1 y) depends on x by computing (3 numerically as a function of fcn/fci2 while keeping k^jkvi 
fixed. We find that h(x,y) depends linearly on x when x^> 1 and goes to zero at x — 1 as (x ~ l) 3 / 2 . We choose the 
fitting function to have the form 



h{x,y) = b(y) 



(x - if 2 
y + (x-iy/ 2 



(25) 



which satisfies these asymptotic requirements. 

The next step is to find a fitting function for b(y). This can be accomplished by requiring that (3 for the cubic 
substrate given by Eq. (24) go to the isotropic limit described by Eq. (22) when fci2 = fcn — 2fc44. This requirement 
uniquely specifies that the function b(y) has the form 



b(y) 



p 



y 



i 



(I) 



1/2' 



3/2 +G/-1 + 1) 3 / 2 ' 



(26) 



where the parameter P is the one determined for the isotropic case, na mely P = 9.14 x 10 -11 . 

To check the quality of the approximation given by Eqs. ([m]) and (pfj), we have computed f3 according to this 
approximation and compared it with (3 obtained by numerical calculations for 88 different materials with cubic 
symmetry. The error in the approximated value of (3 never exceeds 6% for this set of materials, and for the vast 
majority of cases it is below 4%. Table [| shows this comparison for several important substrate materials. It should 
be noted that a given error in (3 results in a smaller error in the predicted damping rate, 7, due to the /3-independcnt 
term in Eq. (||). 



VI. CONCLUSION 

In this paper we have applied the general result of our continuum elastic modelS of phonon-mediated adsorbate 
vibrational relaxation to derive analytic formulas for the FT damping rate when the adsorbates form aggregates on 







the surface. We have considered two possible forms of adsorbate aggregation: (i) when the adsorbates are aligned 
along step edges, and (ii) when the adsorbates form well separated islands. We have derived the expression for the 
damping rate of an infinite periodic chain of the adsorbates, and then shown that the damping rate of a finite locally 
periodic chain is described by the same expression for all chains above a small threshold length. For the case of an 
isolated island we have shown that, under mild restrictions on the island geometry, there is a wide range of possible 
FT resonant frequencies for which the damping rate is essentially the same as for a periodic overlayer with coverage 
equal to the island's local coverage. 

Finally, to facilitate the use of our previously derived expression for the FT damping rate at low coverage, we have 
developed an analytical formula which approximates the parameter /3 appearing in the damping law. The formula 
is valid for both isotropic and cubic materials and achieves remarkable accuracy. We have tested the formula on 
88 different cubic materials (all of whose elastic constants were readily available) and have found that our analytic 
approximation predicts (3 with accuracy better than 6%. 
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